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Abstract. The F5 subgiant Procyon A (a CMi, HR 2943) was observed with the Coralie fiber-fed echelle spectro- 
graph on the 1.2-m Swiss telescope at La Silla in February 1999. The resulting 908 high-accuracy radial velocities 
exhibit a mean noise level in the amplitude spectrum of 0.11 ms~^ at high frequency. These measurements show 
significant excess in the power spectrum between 0.6-1.6 mHz with 0.60 ms~^ peak amplitude. An average large 
spacing of 55.5 fiRz has been determined and twenty-three individual frequencies have been identified. 
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1. Introduction 

The success of helioseismology encourages corresponding 
investigations on other stars. The measurement of frequen- 
cies and amplitudes of p-mode oscillations for solar-like 
stars provides an insight into the internal structure and 
is currently the most powerful constraint on the theory of 
stellar evolution (see Guenthcr & Demaraue ll993|l . 

Solar-like oscillation modes generate periodic motions 
of the stellar surface with periods in the range of 3-30 
minutes but with extremely small amplitudes. Essentially, 
two methods exist to detect such a motion: photome- 
try and Doppler spectroscopy. In photometry, the oscil- 
lation amplitudes of solar-like stars are within 2-30 ppm, 
while they are in the range of 10-150 cms~^ in radial 
velocity measurements. Photometric measurements made 
from the ground arc strongly limited by scintillation noise. 
To reach the needed accuracy requires observations made 
from space. In contrast, Doppler ground-based measure- 
ments have recently shown their ability to detect oscilla- 
tion modes in solar-like stars. Since 1999, solar-like os- 
cillations have been detected in a growing list of main 
sequence and subgiant stars (see the review by Bouchy & 
Carrier 12003) • 

Procyon is the first northern hemisphere candidate for 
the search for p-mode oscillations. Several groups have 
made thorough attempts to detect the signature of oscil- 
lation modes on this bright F5 subgiant star. The first 
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convincing result was obtained by Martic et al 
with observations conducted with the Elodie spectro 
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* Based on observations collected with the Coralie echelle 
spectrograph on the 1.2-m Euler Swiss telescope at La Silla 
Observatory, ESO Chile. 



graph. These Doppler observations led to the detection 
of p-modes in the 0.6-1.6 mHz frequency range with am- 
plitudes 2 times solar and a most likely frequency spac- 
ing equal to 55 /iHz, but no individual frequency was 
identified. On the theoretical side, several stellar evolu- 
tion models of Procyon have been computed predicting a 
large spacing between 50 and 60 /iHz (Barban et al. 119991 
Chaboyer et al. 1999 , Di Mauro & Christensen-Dalsgaard 
nUUU Provost et al. 1^001 and Kervella et al.^UU^. 

In this paper, we report Doppler observations of 
Procyon A made with the Coralie spectrograph, well 
known for the characterization of p-modes on a Cen A 
(Bouchy & Carrier [2001ll200^ . These new measurements, 
made independently with a different calibration method 
than the one used by Martic et al. H1999(l on Elodie (si- 
multaneous thorium method instead of a tunable Fabry- 
Perot illuminated by a white source), confirm the detec- 
tion of p-modes and enable the identification of twenty- 
three individual mode frequencies. 

The data reduction and the acoustic spectrum analysis 
of Procyon A are presented in Sect. 2 and 3, respectively. 
The conclusion is given in Sect. 4. 

2. Observations and data reduction 

Procyon was observed over a campaign of eleven nights 
(8-18 February 1999) with Coralie, the high-resolution 
fiber-fed echelle spectrograph mounted on the 1.2-m 
Swiss telescope at La SiUa (ESO, Chile). 

A description of the spectrograph and the data reduc- 
tion process is presented in Bouchy & Carrier H2001|l and 
Bouchy et al. fjSOOli . The exposure times were 60s with 
110 s dead-times in between. In total, 908 optical spec- 
tra were obtained with a signal-to-noise (S/N) ratio in 
the range of 200 to 300 at 550 nm. For each night, radial 
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Fig. 1. Radial velocities of Procyon relative to a reference 
spectrum taken during each night (the best spectrum of 
the night). The dispersion (which includes the noise and 
the oscillations) reaches 3.1 ms^^. 

Table 1. Distribution and dispersion of Doppler measure- 
ments. 



Date 


Nb spectra 


Nb hours 


a[m s ''] 


1999/02/08 


76 


6.22 


2.83 


1999/02/09 


89 


5.43 


2.92 


1999/02/10 


98 


6.36 


2.97 


1999/02/11 


93 


6.30 


2.95 


1999/02/12 


97 


4.00 


3.45 


1999/02/13 


95 


6.04 


3.03 


1999/02/14 


79 


5.78 


2.59 


1999/02/15 


68 


5.72 


3.17 


1999/02/16 


68 


5.84 


2.56 


1999/02/17 


75 


5.65 


3.32 


1999/02/18 


70 


5.19 


2.61 



velocities were computed relative to the highest signal-to- 
noise ratio optical spectrum obtained in the middle of the 
night (when the target had the highest elevation) . The ra- 
dial velocity measurements are shown in Fig. ^ and their 
distribution and dispersion are listed in Tabled The dis- 
persion of these measurements reaches 3.1 ms^^. 



3. Stellar power spectra analysis 

In order to compute the power spectrum of the veloc- 
ity time series, we use the Lomb-Scargle modified algo- 
rithm (LombUnZni ScargleHnHS giving directly a power 
in m^ s~^ without the need for correction by other nor- 
malization factors. The time scale gives a resolution of 
1.13 /uHz. Fig. |2] shows this periodogram. Typically for 
such a power spectrum, the noise has two components: 

— At high frequencies it is flat, indicative of the Poisson 
statistics of photon noise. The mean white noise level 
CTpow calculated between 2 and 2.5 mHz is 0.015 m^ s^^, 
namely 0.11 ms~^ in amplitude. With 908 measure- 
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Fig. 2. Power spectrum of the radial velocity measure- 
ments of Procyon with the observational window response. 
The white line indicates the noise in the power spectrum. 

ments, this high frequency noise corresponds to aav = 
= 1.85 ms^i. 

— Towards the lowest frequencies, the power should scale 
inversely with frequency squared as expected for in- 
strumental instabilities. However, the computation of 
the radial velocities introduces a high pass filter. 
Indeed, the radial velocities were computed relative 
to one reference for each night and the average radial 
velocities of the night fixed to zero (see Sect. 2). This 
results in an attenuation of the very low frequencies 
which can be seen on Fig. [3 

The power spectrum presents an excess in the range 0.6- 
1.6 mHz. The combined noise has a value decreasing from 
0.027 to 0.016 m^ s~^ in the above-mentioned interval (see 
Fig. 121) . The noise has been determined by fitting a func- 
tion of the type 1 /frequency^ without considering the at- 
tenuated values at very low frequencies. Note that the fil- 
tering induced by the radial velocities computation does 
not infiuence the frequency of the peaks in the range 0.6- 
1.6 mHz, but could slightly change their amplitudes. The 
amplitude of the strongest peaks reaches 0.60 ms~^, corre- 
sponding to a signal to noise of 5 (in the amplitude spec- 
trum). This amplitude is estimated as the height of the 
peaks in the power spectrum with a quadratic subtrac- 
tion of the mean noise level. To investigate if the peaks 
are due to noise or are p-modes, we have conducted sim- 
ulations in which we analyzed noise spectra containing no 
signal. For this purpose, a velocity time series is build, us- 
ing the observational time sampling and radial velocities 
randomly drawn by assuming a Gaussian noise (Monte- 
Carlo simulations). The simulated noise does not need to 
include a spectral dependance, since the threshold deter- 
mination is made relative to a (see below). Indeed, the 
spectral dependance is already included in the observed 
a (noise), which varies with frequency in the power spec- 
trum (see Fig.|2Jl. The amplitude spectrum of this series is 
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then calculated and peaks with amplitude greater than 3, 
4 and 5 a are counted; note that a peak and its aliases are 
only counted once. The whole procedure is repeated 1000 
times to ensure the stability of the results. In this way, 
we find that the number of peaks due to noise with an 
amplitude larger than 3 a varies between and 12 in the 
range 0.6-1.6 mHz, with a mean value of 2.6 and a stan- 
dard deviation of 3.4. For 4 cr, the number of peaks due 
to noise varies between and 2 in the range 0.6-1.6 mHz, 
with a mean value of 0.0 and a standard deviation of 0.3. 
No peaks due to noise are expected with an amplitude 
larger than 5 a. These very convincing results show that 
the power excess is due to p-modes. 

3.1. Search for a comb-like pattern 

In solar-like stars, p-mode oscillations are expected to 
produce a characteristic comb-like structure in the power 
spectrum with mode frequencies Vn,i reasonably well ap- 
proximated by the simplified asymptotic relation (Tassoul 

mm- 
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Fig. 3. Autocorrelation of the power spectrum of Procyon 
with a threshold of 0.25 m^s~^. The dashed lines corre- 
spond to 25.5 /^Hz and its aliases, while the dotted ones 
indicate the value of 31 /iHz and its aliases (see Sect. l3.l|) . 
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(1) 



with Ai^ = {vn.i — i^n-i.e) and 

The two quantum numbers n and I correspond to the 
radial order and the angular degree of the modes, respec- 
tively. Ai/ and bv^^i are the large and small spacing. For 
stars whose disk is not resolved, only the lowest-degree 
modes (I < 3) can be detected. In the case of stellar ro- 
tation, the degeneracy of the modes is lifted and p-modes 
need to be characterized by a third quantum number m 
called the azimuthal order: 



(2) 



with —£ < m < £ and flnim the averaged angular veloc- 
ity over the whole volume of the star. 

One technique, commonly used to search for period- 
icity in the power spectrum, is to compute its autocor- 
relation. To reduce the uncertainties due to the noise, 
only peaks greater than 0.25 m^ s~^ in the power spectrum 
(corresponding to a S/N of 4 in the amplitude spectrum at 
1.6 mHz) in the frequency range 0.6-1.6 mHz have been 
used to compute the autocorrelation shown in Fig. |3| 

The two strongest peaks at 11.5 and 23 /iHz corre- 
spond to the daily alias. Two other peaks at 25.5 and 
31 /LtHz and their daily aliases can be identified. According 
to stellar evolution models of Procyon (Barban et al. 119991 
Chaboyer et al. 119991 Di Mauro & Christensen-Dalsgaard 
EUniland Provost et al. lMT^ . the large spacing Av is ex- 
pected to lie between 50 and 60 /iHz. Fig. 3 exhibits no 
peak having a significant amplitude in this frequency in- 
terval. It is not surprising that a clear peak corresponding 
to the value of the averaged large spacing Ai/ cannot be 
identified with an autocorrelation, since theoretical mod- 
els predict that Ai/ varies sensitively (several /iHz) with 



the frequency between 0.6 and 1.6 mHz. Thus, if only a 
few consecutive oscillation modes are observed, the dif- 
ferent values of the large spacing will not contribute to 
raising the amplitude of a single peak in the autocorre- 
lation, but will instead give a succession of small peaks 
around its averaged value. 

In the case of Procyon, it is much more interesting 
to look for characteristic spacings between modes with 
different angular degrees. Indeed, theoretical models pre- 
dict that, contrary to the large spacing Aiy, these spacings 
remain approximately constant with the frequency. The 
peak at 25.5 /iHz can be identified as the spacing between 
the modes £ = 1 and £ = with the same radial order 
and the spacing between £ = 2 and £ = 1, also with the 
same radial order. The peak at 31 /iHz results from the 
spacing between the modes {£ = 0, n = hq) and {£ = 1, 
n — no — 1), as well as the spacing between {£ = 1, n ~ uq) 
and {£ ~ 2, n = uq — 1). The autocorrelation spectrum 
thus suggests that the large spacing must be close to the 
sum of these two peaks, i.e. close to 56 /^Hz. 

3.2. Echelle diagram 

To identify the individual frequencies of the modes, we 
made echelle diagrams corresponding to values of the av- 
eraged large spacing around 56 /iHz. 

The frequencies were extracted using an iterative al- 
gorithm that identifies the highest peak between 0.6 and 
1.6 /iHz and subtracts it from the time series. Note that 
because of the stochastic nature of solar-like oscillations, 
a timestring of radial velocities cannot be expected to 
be perfectly reproduced by a sum of sinusoidal terms. 
Therefore, using an iterative clean algorithm to extract 
the frequencies can add additional peaks with small am- 
plitudes due to the finite lifetimes of the modes that we 



4 P. Eggenberger et al.: Solar-like oscillations in Procyon A 

0.6 I 1 ^I^ T— I 1 J "-^^ 1^ TT , , , , ^1 n| 1 , , ^1 1 n , rn f 




600 BOO 1000 1200 1400 1600 

Frequency [mHz] 

Fig. 5. Power spectrum of Procyon with the twenty-eight extracted frequencies indicated by dotted hnes. The identi- 
fication of each extracted frequency is given in Table |2 
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Fig. 4. Top: Original power spectrum for the eleven 
nights of observation. Bottom: Cleaned power spectrum: 
all peaks listed in Table[21have been removed. The contin- 
uous, dotted and dashed lines indicate an amplitude of 5 ct, 
4ct and 3ct, respectively. Numerous peaks are still present 
below 3 (T, since no peaks have been cleaned below this 
threshold. These peaks can be due to p-mode oscillations 
and noise or have been artificially added by the extraction 
algorithm due to the finite lifetimes of the modes 

do not know. Nevertheless, the iterative algorithm ensures 
that one peak and its aliases with an amplitude above a 
given threshold is only extracted once. To avoid extract- 
ing artificial peaks with small amplitudes added by the 
iterative algorithm, the choice of this threshold is impor- 
tant. In the case of Procyon, we iterated the process until 
all peaks with an amplitude higher than 3 ct in the am- 
plitude spectrum were removed (see Fig. QJ. Note that 



Table 2. Identification of extracted frequencies. Some fre- 
quencies can be either split £ — 2 mode or due to noise. 
The number of peaks due to noise is in agreement with the 
simulations described in Sect. 13 which predict 2.6 ± 3.4 
noise peaks. 



Frequency 




Mode ID 


S/N 


651.5 - 11.57 = 


= 639.9 


£ = 2 


3.2 


630.8 + 11.57 = 


= 642.4 


£ = 


3.9 




662.7 


noise 


3.2 


683.5 + 11.57 = 


= 695.1 


£ = 


4.7 




720.6 


£= 1 


4.5 




797.9 


£ = 2 


4.1 




799.7 


£ = 2 


4.8 


791.8 + 11.57 = 


= 803.4 


£ = 


3.2 




828.5 


£= 1 


3.0 




835.4 


noise 


3.6 




856.2 


£ = 2 


4.4 




859.8 


£ = 


4.3 




911.4 


£ = 2 


3.2 


929.2 - 11.57 = 


= 917.6 


£ = 


3.9 


929.2 + 11.57 = 


= 940.8 


£= 1 


3.9 


1009.7 - 11.57 = 


= 998.1 


£= 1 


4.5 




1027.1 


£ = 


3.4 


1123.3 + 11.57 = 


1134.9 


£ = 2 


4.8 




1137.0 


£ = 2 


4.0 


1131.1 + 11.57 = 


1142.7 


£ = 


3.0 




1192.4 


£ = 2 


3.0 


1186.0 + 11.57 = 


1197.6 


£ = 


3.4 


1234.8 + 11.57 = 


1246.4 


£ = 2 


4.9 




1251.8 


£ = 


3.7 


1265.6 + 11.57 = 


1277.2 


£ = 1 


5.0 




1337.2 


noise 


4.0 




1439.0 


noise 


3.3 




1559.5 


£ = n 


3.7 



peaks with amplitudes below the 3 a threshold were not 
considered not only to avoid extracting artificial peaks 
added by the iterative algorithm, but also because they 
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Table 3. Mode frequencies (in fiHz) of Procyon. 



£ ^ 


1 = 1 


£ ^ 2 






639.9 


642.4 






695.1 


720.6 


797.9 / 799.7 


803.4 


828.5 


856.2 


859.8 




911.4 


917.6 


940.8 
998.1 




1027.1 




1134.9 / 1137.0 


1142.7 




1192.4 


1197.6 




1246.4 


1251.8 


1277.2 





were too strongly influenced by noise and interactions be- 
tween noise and daily aliases. This threshold, which en- 
sures that the selected peaks have only a small chance 
to be due to noise, gave a total of twenty-eight frequen- 
cies (see Table El). Because of the daily alias of 11.57 /iHz 
introduced by the monosite observations (see the obser- 
vational window response in Fig. |2Jl, we cannot know a 
priori whether the frequency selected by the algorithm is 
the right one or an alias. Thus, we considered that the fre- 
quencies could be shifted by ±11.57 /iHz. We then made 
echelle diagrams for different large spacings between 50 
and 60 /iHz until each frequency could be identified as 
an^ = 0, ^ = 1, £ ^ 2 mode or attributed to noise (see 
Table 121). In this way, we found an averaged large spacing 
of 55.5 /iHz. 

At 799 and 1136 /iHz, cither the £ = 2 modes are 
split by the rotation or the second peaks are due to 
the noise. AUende Prieto et al. H2U(J2[) recently estimated 
vs'mi — 3. 16 ±0.50 km s~^ but mentioned that the correct 
value is probably close to 2.7kms~^, given that the value 
of 3.16kms~^ may be slightly overestimated due to the 
finite numerical resolution of their convection simulation. 
Using vs'mi = 2.7kms~^ and i = 31.1±0.6° (Girard et al. 
I2000|) . the rotational splitting expected for Procyon A is 
about 0.6 /iHz. For £ — 2 modes, this splitting can thus re- 
sults in differences between two modes as high as 2.4 /iHz, 
which is the frequency separation between m — ^2 and 
ni — +2 modes. Of course, the quality of our data is not 
good enough to unambiguously determine this rotational 
splitting. The possible £ = 1 peak at 1559.5 /iHz is re- 
jected, as it is far from the other modes with a doubtful 
identification. Moreover, the frequency of 1027.1 /j,Hz is 
difficult to identify; we identify it as a £ = mode, but it 
could also he a, £ — 2 mode. 

The echelle diagram showing the twenty-three identi- 
fied modes is shown in Fig. El The frequencies of the modes 
are given in Table 
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Fig. 6. Echelle diagram for the frequencies listed in 
Table O Dots, triangles and squares correspond respec- 
tively to modes identified as ^ = 0, £ = 1 and £ = 2. The 
size of the symbols is proportional to the signal-to-noise 
ratio of the mode. 

3.3. Amplitudes of the modes 

Concerning the amplitudes of the modes, theoretical com- 
putations predict oscillation amplitudes for a 1.5 Mq star 
like Procyon A between 1 and 2 ms^"'^, with mode life- 
times of the order of a few days (Houdek et al. ll999f) . The 
observed amplitudes of 0.6 ms""'^ are then lower than ex- 
pected. This disagreement can be partly explained by the 
lifetimes of the modes. Indeed, the oscillation modes have 
finite lifetimes, because they are continuously damped. 
Thus, if the star is observed for a time longer than the 
lifetimes of the modes, the signal is weakened due to the 
damping of the modes and to their re-excitation with a 
random phase. 

4. Conclusion 

The radial velocity measurements of Procyon A, obtained 
over 11 nights, show a significant excess in the power spec- 
trum between 0.6-1.6 mHz, centered around 0.8 mHz, 
with a peak amplitude of 0.60 ms~^, revealing solar-like 
oscillations with a large spacing of 55.5 /iHz. Our results 
confirm the values of the large spacing and of the ampli- 
tudes of the modes found by Martic et al. H1999|) . 

Moreover, we presented the identification of twenty- 
three individual frequencies. Note that this identification 
is seriously complicated by the presence of daily aliases. 
Interactions between real peaks and aliases can slightly 
shift the frequencies of the modes. This is particularly 
true for the £ — and £ — 1 modes which are separated by 
25.5 /iHz: the interactions between the second alias of the 
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radial modes (separated from the real peak by 23.14 fjHz) 
and the £ — 1 modes can slightly shift the observed £ — 
1 frequencies, although the different peaks are resolved. 
Inversely, the second alias of the £ = 1 modes can influence 
the observed frequencies of the radial modes. Nevertheless, 
we think that this identification is reliable, since it is the 
only one that can explain all the highest peaks in the 
power spectrum. 

To obtain more accurate parameters and to thoroughly 
test the physics of the models, it is not only important to 
obtain better data, but also to eliminate the mode iden- 
tification ambiguity due to the aliases. The space mission 
Most fMatthews i2000.) launched on 30 June 2003 wiU ob- 
serve Procyon A during a whole month and should there- 
fore provide severe constraints on this star. 
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